Cosmic strings from preheating 
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We investigate nontliermal phase transitions that may occur after post-inflationary preheating in 
a simple model of a two-component scalar field with the effective potential X{4>i — v^)^/4, where 
01 is identified with the inflaton field. We use three-dimensional lattice simulations to investigate 
the full nonlinear dynamics of the model. Fluctuations of the fields generated during and after 
preheating temporarily make the effective potential convex in the (/ii direction. The subsequent 
nonthermal phase transition with symmetry breaking leads to formation of cosmic strings even for 
V S> 10^® GeV. This mechanism of string formation, in a modulated (by the oscillating field 01 ) 
phase transition, is different from the usual Kibble mechanism. 
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The theory of cosmological phase transitions proposed 
in 1972 by David Kirzhnits gradually became one of the 
most essential parts of modern cosmology Origi- 
nally it was assumed that such phase transitions occur 
in a state of thermal equilibrium when the temperature 
decreased in the expanding universe. Recently it was 
found that large fluctuations of scalar and vector fields 
produced during preheating after inflation ||] may lead 
to specific nonthermal phase transitions which occur far 
away from the state of thermal equilibrium j^. To in- 
vestigate these phase transitions one should study the 
self-consistent nonlinear dynamics of quantum fluctua- 
tions amplifled by parametric resonance. This is a very 
complicated task. Fortunately, fluctuations of Bose fields 
generated during preheating have very large occupation 
numbers and can be considered as interacting classical 
waves, which allows one to perform a full study of all 
nonlinear effects during and after preheating using lattice 
calculations B. These calculations, as well as analytical 
estimates have already shown that fiuctuations can 
grow large enough for cosmologically interesting phase 
transitions to occur. Lattice calculations can be used to 
directly simulate nonthermal phase transitions and for- 
mation of topological defects. This made it possible to 
go beyond the early attempts to study such phase tran- 
sitions numerically |^ , which neglected crucial backreac- 
tion effects beyond the Hartree approximation. 

We made a series of lattice simulations of nonthermal 
phase transitions, focusing primarily on the possibility of 
generation of topological defects ||l^,|lll . The model with 



one scalar field with the potential \{(jr 



74 has a 



broken discrete symmetry (f> ^ —(p. In this model we 
have found formation of domain structure pc| ]. If the 
inflaton field strongly couples to another scalar field x, 



the phase transition is strongly first order; we observed 
formation of bubbles of the new phase |jll[] . 

In this Letter we will study string formation in the the- 
ory of a two-component scalar field with the effective 
potential A((/)^ — v^)^/4, where 0^ = "12^1 ■ This model 
has 0(2) rotational symmetry and allows for string for- 
mation |l^ (for early conference reports see [0). Re- 
cently Kasuya and Kawasaki reported formation of de- 
fects in 2d lattice simulations for v < 3 x 10^^ GeV @, 
which confirmed our general conclusion concerning the 
phase transition in this model ||l^,|l^. However, scat- 
tering of particles, as well as the nature of topological 
defects in two dimensions and in three dimensions are 
quite different. In particular, there are no strings in this 
model in two dimensions. We have found that in the 
realistic case of three dimensions the generation of fluc- 
tuations is much more effective, and string production is 
possible even if v is as large as 6 x 10^^ GeV. 

In the model with the two-component field (pi we can 
always rotate fields in such a way that initially 02 — 0, 
and the field (pi plays the role of the classical oscillating 
inflaton field. We denote by 0(0) the value of the field 
01 at the moment when inflation ends and the inflaton 
field begins to oscillate. It is convenient to work with the 
rescaled conformal time r, where ^f\(p{{))dt = a{T)dT, 
a(0) = 1, and perform conformal transformation of the 
fields, ip — 0(T)a(T)/0(O). We also wfll use the rescaled 
spatial coordinates x —^ \/A0(O)x. The initial conditions 
at the beginning of preheating are determined by the 
preceding stage of inflation. We deflne the beginning of 
preheating as the moment r = when the velocity of 
the field in conformal time is zero. This happens when 
0(0) « 0.35Mpi and a(r) w 0.51r -fl g. 

In the new variables, the equations of motion become 
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(pi - v'^lpi + {ip^ 
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where v = v/(/)(0). Equation (|l]) contains only one pa- 
rameter, V. The coupUng constant A is hidden in the 
initial conditions for fluctuations; these were chosen as 
described in Ref. Q|. 

The full nonlinear equations of motion were solved 
numerically directly in the configuration space. The com- 
putations were done on 128'^ lattices with the box size 
L = IGtt, with the expansion of the universe assumed 
to be radiation dominated. There are several important 
quantities that can be measured in the simulations dur- 
ing the evolution of the system. One can define the zero 
mode ipai — {(pi {x, r)) for each of the components i = 1, 2 
of the field ip, and variances {(Sipi)'^), which measure the 
average magnitude of fluctuations Sipi = ipi{x,T) — ipQi. 
Since the system is homogeneous on large scales, averages 
in these relations can be understood as volume averages. 
This is equivalent to taking an average over realizations 
of the initial data (the ensemble average). The averaged 
quantities depend only upon time and do not depend 
upon spatial coordinates x. 

We performed calculations for various v and A. In 
the beginning we present results for v = 3 x 10^^ GeV 
and A = lO^^'^. In the last part of the paper we also 
present result for different values of v and for A = 10~^^. 
The behavior of the zero modes (jjoi = ^nd of the 

variances, {{(f>i — <j>oi)'^), for each of the field components 
is shown in Fig. |^. 
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FIG. 1. Variances of the field components 02 and <jf>i (up- 
per and lower solid curves respectively) as well as the zero 
modes of these fields (upper dotted curve is {(j>i) and lower 
dotted curve is {(j)2)) are shown as functions of conformal time. 

Let us explain the evolution of the fluctuations. For 
such a small value of v, preheating is completed at 
r « 100, well before the phase transition with symmetry 
breaking, which occurs at r « 240. During the stage of 
parametric resonance r < 100, the term u^a^ in (|l]) can 
be neglected. The equations for the mode functions of 
fluctuations in the directions of 0i and 02 are 



where the background inflaton oscillations (Poi{t) are 
given by an elliptic function. The resonance parameter 
is different for different components: qi = 3 and q2 = 1. 
For both components there is a single instability band, 
but the locations of the bands and the strengths of the 
resonance are different jlj]. The resonance in the "in- 
flaton" direction ipi is weak, the maximal value of the 
characteristic exponent of the fluctuations ipik oc e'^^'^ is 
/ii « 0.036; the resonance in the perpendicular direction 
ip2 is much stronger and broader, /Lt2 ~ 0.147. 

However, the actual growth rate of the fluctuations tpik 
is not smaller but larger than that of the ip2k fluctuations. 
Indeed, the cross-interaction term Lp1ip2 leads to the pro- 
duction of Stpi fluctuations in the process of rescattering 
Q. As a result, the amplitude ipi^k grows as ipik oc e^'^^'^, 
where 2/i2 = 0.294 |12], as clearly seen in Fig. 0. This 
growth is much faster than in the 2d lattice simulations 
of Ref. |l|. 

We see that during the time interval between comple- 
tion of preheating and the phase transition, 100 ^ t < 
240, the zero mode of (f>i decreases faster than the vari- 
ances (both are decreasing due to the expansion of the 
universe, but in addition the zero mode continues to de- 
cay into field fiuctuations) , and by the time r « 240 the 
field variance is as big as its zero mode. 
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FIG. 2. Time dependence of the zero mode of (j>\. 

Time interval 100 ^ 240 is dominated by rescat- 
tering of all modes [Q. During this time interval, the 
symmetry (at least in one direction) is restored by large 
fiuctuations. To see this, we study the time dependence 
of the zero mode of 0i, which is shown in Fig. ||. We 
immediately see some peculiarity. In the potential of the 
form {(fP' — v^)^ the field cannot oscillate near (f) — with 
amplitude smaller than 0c = v^v. This critical value 
of the amplitude is shown in Fig. ^ by the dashed line, 
and we see that the field does oscillate with an amplitude 
smaller than (pc- At some point, the amplitude of the os- 
cillations even becomes smaller than v, i.e. the field is 
oscillating on the top of the tree potential without rolling 
down to its minimum. This means that the tree poten- 
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tial is significantly altered by the interaction with the 
background of created fluctuations, and the symmetry is 
restored. 

We can reconstruct the effective potential Vcg{(p) using 
the already calculated function (poi^r) — (ipi) and the 
definition 

ipoi + dVcs/dcpoi = . (3) 

This definition of the effective potential is perfectly legit- 
imate in the case when corrections to terms with the time 
derivative are small, which appears to be the case. (Ex- 
pansion of the universe will be properly included when 
we work in the conformal coordinates, as in Eq. (j^), and 
then "rotate" potential back to the synchronous frame 
and to the original fields.) This method allows one to 
find VcS up to a constant. The reconstruction along ipi 
direction is shown for several moments of time in Fig. |^. 
At T w 90 the effective potential, shown by dots, coin- 
cides to a very good accuracy with the tree potential, 
shown by the solid line. This moment is not far away 
from the end of the exponential growth of fluctuations, 
which occurs at r « 100, but fluctuations at r « 90 are 
still rather small. At r « 110 fluctuations are large, and 
the effective potential is completely different: it has only 
one minimum. Note that the potential is slightly asym- 
metric with respect to the change of sign of the field 0i. 
This happens because in our simulations we are sampling 
the effective potential in the time-dependent background, 
when the points with positive and negative (jj correspond 
to different moments of time. The "instantaneous" po- 
tential would be exactly symmetric, with a minimum at 
01 = 0, and symmetry is restored, at least in the 0i 
direction. 
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FIG. 3. Reconstruction of the effective potential in (f>i 
direction at several moments of time. Dots correspond to 
r r; 90, diagonal crosses to t « 110, larger stars to r ~ 180 
and smaller stars to r ~ 215. Solid line is the tree potential. 

Still, the existence of an oscillating zero mode 0i im- 
plies that symmetry between states with 0i > and 
01 < is not exact. Thus one wonder whether one can 
say that symmetry is restored if the effective potential 



has the minimum at 0i = 0, or one should not say so un- 
til the amplitude of oscillations of the field 0i completely 
vanishes. 

In our opinion, if the effective potential has a minimum 
at 01 = 0, this implies that there is no spontaneous sym- 
metry breaking, so in this sense the symmetry 0i ^ — 0i 
is restored. From a more pragmatic point of view, the 
issue of symmetry restoration and subsequent symme- 
try breaking is important mainly because it is related 
to production of topological defects. Therefore instead 
of debating which definition of symmetry restoration is 
better, we will try to find whether the topological defects 
are produced. Indeed we will see that cosmic strings are 
produced as a consequence of preheating. 

At T « 240 the field fluctuations become diluted by the 
expansion of the universe to the extent sufficient for sym- 
metry breaking. This event can be seen both in Fig. |l| 
and in Fig. ^. At that time, the zero mode nearly van- 
ishes. This happens not because of symmetry restoration, 
but because the field rolls down to the 0{2) symmetric 
valley of minima of the effective potential in all possible 
directions in the field space, which, after averaging, gives 

(0l) < V. 




FIG. 4. The joint probability distribution function of (j>i 
and (j}2 before and after the phase transition. Fields (pi and 
(l>2 are shown in units of v. 

A useful quantity to consider is the probability distri- 
bution function P(0i, 02, t) shown in Fig. ^. The first of 
these two figures shows that at r ~ 200 the maximum of 
the probability distribution oscillates near 0i = in the 
01 direction. This distribution has only one maximum 
in the 0i direction. When this maximum approaches 

01 = 0, the distribution P(0i,02,t) is approximately 
symmetric with respect to the change 0i — > — 0i. On 
the other hand, there are two maxima of the probability 
distribution with respect to the field 02. They are con- 
centrated near 02 ~ ±v, which means that the symmetry 

02 —02 is broken. 

This implies that the universe at that time becomes di- 
vided into domains filled with the field 02 ~ ±v. These 
domains are separated by two-dimensional domain walls, 
the surfaces where 02 = 0. When the distribution of 
the scalar field 0i oscillates, the center of this distribu- 
tion moves, but if it is wide enough, there always will 
be two-dimensional surfaces where 0i = 0. Intersection 



3 



of these surfaces with the domain walls 02 = form 
strings, on which 0i = 02 = 0. One can easily find out 
that when one moves around such a string, the phase 
a = arccos(0i/|0|) changes by 2n. Thus these strings 
are topologically stable. Most of them form string loops 
which move (expand, shrink, and expand again) when the 
distribution of the field 0i oscillates. This mechanism of 
string formation is different from the Kibble mechanism. 

Gradually the amplitude of fiuctuations of the fields 0i 
decreases, and symmetry 0i — s- —<j)i also breaks down. 
The choice of direction of the symmetry breaking will 
depend on the shape of the effective potential, but also 
on the amplitude and phase of the oscillations of the zero 
mode (0i), and on the width of the probability distribu- 
tion P{(j)i, <j>2,t). If the width of the probability distribu- 
tion P{(j)i, (j)2,t) in the 4>i direction is sufficiently large, 
then this phase transition modulated by the oscillations 
of the field 0i may lead to formation of long strings, 
which may have interesting cosmological consequences. 

After the modulated phase transition, which occurs at 
T — 240, we see an 0{2) symmetric ring in the probabil- 
ity distribution P{(j)i, (j)2,t), superimposed with a peak 
in some random direction (which does not coincide ei- 
ther with 01 or with 02), see the second figure in Fig. ^. 
The existence of the ring shows that the absolute value 
of the field after spontaneous symmetry breaking is close 
to V, and that in different points of space the vector 
(01,02) looks in all possible directions, which indicates 
the presence of strings. The peak along the random direc- 
tion represents additional spontaneous symmetry break- 
ing, which appears because it is energetically preferable 
for all vectors (0i, 02) to look in the same direction. With 
time, the ring will disappear, and the width of the peak 
will decrease. 




Strings can be detected directly. For that purpose, we 
had plotted 3D coordinates of the grid points where |0| is 
close to zero. A series of such images describing several 
different stages of this process is shown in Fig. |5[ Strings 
are clearly seen, and we can observe the formation of one 
large loop. 

At times 100 < r < 240, while the zero mode was 
still oscillating, we had observed abundant formation and 
subsequent annihilation of string loops. Many strings 
appear when the oscillating zero mode (0i(i)) passes 
through 01 = 0. Then they disappear when (0i) grows 
and appear when it becomes small again. Fig. |^, on the 
other hand, shows the behavior of strings after the phase 
transition, when the field 0i is not capable of rolling over 
the top of the effective potential. At this stage the new 
strings are not produced, and the old ones move slowly, 
changing their position because of the string tension. 

It is interesting that at time r = 280 there is only 
one string loop in the integration box (upon account of 
the periodic boundary conditions), see Fig. ||. (A single 
spatial configuration, like this one, of course corresponds 
to a particular realization of random initial conditions.) 
This loop stays intact, up to small vibrations, for a long 
time, being in quasi-equilibrium. However, later on, dif- 
ferent segments of the loop quickly approach each other 
and reconnect at r = 328, forming another loop config- 
uration at T — 330. This final loop collapses almost to 
a point, bounces once, collapses again and disappears.^ 
This process of reconnection and final collapse confirms 
that loops are dynamical strings, not accidental lines in 
space. 

To find the range of v for which strings are produced 
we had run simulations with different values of v in the 
interval from 4 x 10^^ to 10^^ GeV. In all these cases 
we observed formation of string loops. From cosmologi- 
cal perspective, however, the most interesting possibility 
would be formation of an infinite string. One can ex- 
pect that, at least when v is large, the probability of 
formation of long strings is higher in the case when the 
moment of the symmetry breaking nearly coincides with 
the moment when (0i) passes through zero. If that ex- 
pectation is correct, the number of long strings should be 
a non-monotonic function of v. 

To verify this conjecture we plotted the string distri- 
bution at the time At — 10 after the phase transition 
(which happens at different times) for v = 3 x 10^^ GeV, 
5x10^^ GeV, 6 X 10^6 GeV, and 10^^ GeV, for A = 10""^ 
see Fig. |6[ As we see, for v = 5 x 10^^ GeV and 10^'' 
GeV all loops are short, whereas for v = 3 x 10^^ GeV 
and 6 x 10^^ GeV there are many large loops. This in- 
dicates a possibility of formation of a network of infinite 



FIG. 5. The process of string formation for v = 3 x 10^^ 
GeV and A = 10" 



*The movie which shows this process can be see n at 



http:/ /www. physics. purdue.edu/'~-^tkachev/movies. htm] or at 



http://physics.stanford.edu/~linde 
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strings for V = 3 X lO^'^ GeV and 6 x lO^*' GeV. 



mechanism, including infinite ones. We hope to return 
to the discussion of this model in a separate publication. 




FIG. 6. The string distribution at the time Ar = 10 after 
the moment of the phase transition for v — 10^^ GeV, 6 x 10^^ 
GeV, 5 X lO^'' GeV, and 3 x lO" GeV respectively. The plot 
for V = 3 x 10^'' GeV corresponds to a different realization of 
the initial conditions than Fig. 5. 

Additional evidence for (or against) formation of infi- 
nite strings after preheating may be obtained if one uses 
lattices of a greater size. In particular, we performed sim- 
ulations for V = 3 X 10^^ GeV in a larger box {L — 32tt), 
see Fig. 0. At the time shown in this figure, the physical 
size of the box is Lphys ~ (10^)m~^, where ~ 1/VAv 
is the typical thickness of the string. In Fig. we see 
two "infinite" strings and a large string loop. Note that 
on a lattice with periodic boundary conditions "infinite" 
strings can only be created in pairs, because the winding 
number for our initial conditions is zero. 

Also, one may consider models where one may have 
independent reasons to expect production of infinite 
strings. For example, one may add to our model the 
term ^0^x^ describing interaction of the infiaton field 
with the scalar field x with the coupling constant 3> A. 
We have studied this model in Ref. for the case of 
a one-component field and found a first-order phase 
transition. We expect that this result remains valid for 
the two-component field as well, because the main reason 
for the first order phase transition found in [0 was the 
strong interaction of the field (f> with the field x rather 
than the self-interaction of the field (p. If this is indeed 
the case, then nonthermal fluctuations of the field x will 
lead to symmetry restoration with respect to both of the 
fields 4>i and (j>2 ■ These two fields will be captured in the 
local minimum of the effective potential at (/) = until 
the moment of the first order phase transition. In this 
case, there will be strings formed by the usual Kibble 




FIG. 7. The string distribution for v = 3 x 10^*^ GeV in a 
box of a larger size. 

This work was supported in part by DOE grant DE- 
FG02-91ER40681 (Task B) (S.K. and I.T.), NSF grants 
PHY-9219345 (A.L.), PHY-9501458 (S.K. and I.T.), 
AST95-29-225 (L.K. and A.L.), and by the Alfred P. 
Sloan Foundation (S.K.). 



[1] 

[2] 

[3] 

[4] 
[5] 
[6] 



D.A. Kirzhnits, JETP Lett. 15, 529 (1972); D.A. Kirzh- 
nits and A.D. Linde, Phys. Lett. 42B, 471 (1972); Sov. 
Phys. JETP 40, 628 (1974); Ann. Phys. 101, 195 (1976); 
S. Weinberg, Phys. Rev. D9, 3320 (1974); L. Dolan and 
R. Jackiw, Phys. Rev. D9, 3357 (1974). 
L. A. Kofman, A. D. Linde, and A. A. Starobinsky, Phys. 
Rev. Lett. 73, 3195 (1994); L. A. Kofman, A. D. Linde, 
and A. A. Staro binsky, Phys. Rev.D 56, 3258 (1997), 



hep-ph/9704452 



L. Kofman, A. Linde, and A. A. Starobinsky, Phys. Rev. 
Lett. 76, 1011 (1996); I. I. Tkachev, Phys. Lett. B376, 
35 (1996). 

S.Yu. Khie bnikov and I.I. Tk achev, Phys. Rev. Lett. 77, 



219 (1996), |hep-ph/9603378 , 

S. Khiebnikov and I. Tkachev, Phys. Lett. B390, 80 



(1997), |hep-ph/9608458 
S 

1607 



Khiebnikov and I. Tkachev 
(1997 ) 



ph/9708313 



hep-ph/9610477 



Phys. Rev. Lett. 79, 
S. Khiebnikov, |iep 



[7] S. Khi ebnikov and I. T kachev, Phys. Rev. D56, 653 



(1997), iiep-ph/9701423 



5 



[8] 
[9] 

[10] 
[11] 

M. 



T. Pro kopec and T. G . Roos, Phys. Rev. D 55, 3768 
(1997), |hep-ph/9610400 . 

D. Boyanovsky, D. Cormier, H.J. de Vega, R. Holman, A. 
Singh, and M. S rednicki, Phys. Rev. D 56, 3929 (1997), 



hep-ph/9703327; S. Kasuya, M. Kawasaki, Phys. Rev. D 



56, 7597 (1997), [hep- ph/9703354 , 

I. I. Tkachev, S. Khlebnikov, L. Kofman, A. Linde, and 
A. A. Starobinsky, in preparation. 

S. Khlebnikov, L. Kofman, A. Linde, and I. Tkachev, 
"Firs t-order nontherm al phase transition after preheat- 
ing, " |hep-ph/9804425| . 

L. Kofman, a talk at COSM097, Ambleside, 1997, ^ep- 
ph/980222l]; L. Kofman, a talk at the 3rd RESECEU 



Symposium on Particle Cosmology, Tokyo 1997; A.D. 

Linde, talks at COSM097, Ambleside, 1997, and at PAS- 

COS98, Boston, 1998. 
[13] S. Kasuya, M. Kawasaki, "Topological defects formatio n 

after inflation on lattice simulations," liep- ph/980442g| . 
[14] D. Boyanovsky, H.J. de Vega, R. Holman, and J. Salgado, 

Phys. Rev. D 5 4, 7570 (1996); D .L Kaiser, Phys. Rev. D 

56, 706 (1997), |hep-ph/9702244 P. B. Green, L. Kofman, 



A. D. Linde, and A. A. Star obinsky, Phys. Rev. D 56, 
6175 (1997), [hep-ph/9705347. 



6 



